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INTRODUCTION 


This  program  was  developed  at  Thiokol  Corporation  to  assist  in 
the  design  of  elastomeric  bearings. 

It  is  a  simplified  version  of  a  larger  system.  Because  of 
this#  it  has  several  "vestigial  organs"  in  It.  For  example# 
the  material  property  arrays  and  routines  are  fully  set  up  for 
orthotropic  materials  even  though  there  is  no  orthotropic 
capability  in  the  program.  Other  instances  of  this  will  be 
noticed  in  the  program.  These  vestigial  functions  have  some 
advantage  in  that  the  capability  of  the  program  can  be 
expanded  more  easily  if  and  when  desired. 

The  program  has  been  designed  to  be  modular  in  nature.  Thus 
particular  functions  are  restricted  to  certain  modules.  This 
also  aids  in  maintenance  and  modification.  Each  of  the 
modules  is  described  in  later  sections. 

In  general  non-FORTRAN  I/O  has  been  used  in  the  program.  This 
has  two  main  advantages.  First#  it  Is  faster#  In  general  two 
to  four  times  faster#  than  binary  FORTRAN  I/O.  Second#  some 
things  are  possible  in  the  special  routines  that  can  not  be 
done  In  FORTRAN.  For  example#  subroutine  TDRR  uses  the 
NOTE/POINT  capability  to  retrieve  large  records  directly  into 
their  place  in  core  in  an  order  different  from  that  In  which 
they  were  wri tten. 
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module  pescriptiqns 


The  basic  philosophy  of  the  program  has  been  to  try  to 
separate  functions  into  moduiues.  Following  are  descriptions 
of  each  of  the  modules  in  the  program.  These  are  not  intended 
to  be  comprehensive  but/  along  with  the  flow  charts/  to  give 
an  introduction  to  the  system. 


Ili£  C.QNim  SECTION 


This  part  of  the  program  is  the  driver  that  controls  the  logic 
flov/  through  each  of  the  modules. 

The  program  begins  execution  in  a  small  FORTRAN  routine  v/hose 
only  function  is  to  call  the  GTMAIN  routine.  The  GTMAIfl 
routine  uses  a  conditional  getmain  to  get  a  large  block  of 
core  storage  to  be  used  as  a  work  area  by  the  program  and  then 
passes  control  to  subroutine  COMTRL. 

Subroutine  CONTRL  actually  controls  the  logic  of  the  program, 
it  reads  the  title  and  the  general  input  section  and 
calculates  basic  dimensions.  A  set  of  logical  variables 
contained  in  common  CNTRL  specifies  v/hich  modules  will  be 
entered  and  how  the  program  will  loop.  These  logicals  are 
initialized  in  a  block  data  routine  and  then  arc  changed 
according  to  program  options. 


Iil£  INPUT  nmiLL 


This  module  controls  the  reading  of  the  input  data  from  cards. 
All  Input  to  the  program  is  via  subroutine  FREFRM  v/hich  allows 
data  to  be  entered  in  a  free  field  format.  Each  of  the 
routines  in  this  module  calls  FREFRM  and  then  interprets  the 
string  of  data  returned  to  it. 

Subroutine  IIIMOD  Is  the  first  routine  called.  It  sets  up  tiie 
basic  arrays  out  of  the  v/ork  area  and  then  passes  control  to 
subroutine  INPUT. 

Subroutine  INPUT  checks  the  flag  record  against  a  table  of 
valid  flags  and  then  calls  the  proper  routine  to  handle  that 
section  of  input.  At  tlie  end  of  the  input  it  calls  subroutine 
PRGEOM/  v/hich  prints  the  geometry  and  creates  the  geometry  and 
load  files  (IGEOM  and  I  LOAD). 
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This  module  uses  for  scratch  files  two  pairs  of  "fl Ip/flop" 
files,  that  is,  the  data  Is  passed  back  and  forth  betv/een  the 
two  files  with  the  file  names  being  interchanged 
(flip/flopped)  betv/een  each  pass.  One  set  Is  used  for  the 
integer  code  information,  and  one  set  is  used  to  hold  the  load 
Information.  This  technique  avoids  the  use  of  large  In-core 
arrays  to  hold  all  of  this  Information.  It  is  a  trade-off 
between  core  size  and  run  time. 


Ili£  SOLUTION  MODULE 


This  module  creates  the  element  stiffness  matrices  and  force 
vectors,  assembles  them,  and  solves  for  the  displacements. 
The  solution  scheme  used  is  a  Gaussian  scheme,  and  only  the 
semi -bandwi dth  for  a  row  of  elements  is  in  core  at  any  one 
time. 

The  control  routine  for  this  nodule  is  SOLMOD.  It  sets  up  the 
arrays  and  calls  subroutine  FORSOL  for  the  forward  part  of  the 
solution,  and  then  reallocates  the  arrays  and  calls  subroutine 
BAKSUB  for  the  backward  part  of  the  solution,  and  the 
extraction  of  the  displacements.  If  a  large  deformation 
iteration  is  Involved,  then  subroutine  LDCTUL  Is  called  to 
control  that  function. 

Subroutine  FORSOL  loops  through  for  each  element  in  the 
system.  The  geometry  and  load  files  are  read  for  each  row  of 
nodes,  and  the  loads  are  adjusted  If  there  Is  an  Incremental 
loading  situation.  Material  properties  are  calculatetl  for 
each  element  and  Interpolated  for  If  there  Is  strain  or 
thermal  dependence.  A  logical  array,  LSTP,  Is  set  up  v/hlch 
contains  a  reference  of  flags  giving  the  degrees  of  freedom  on 
each  node;  this  array  Is  used  In  the  assembly  of  the  solution 
matrix  and  the  extraction  of  the  roots.  The  displacement 
boundary  conditions  and  nodal  point  forces  are  applied  to  each 
element  stiffness  matrix.  The  solution  matrix  is  assembled 
and  partially  reduced  for  a  row  of  elements  at  a  time  and  the 
equations  arc  written  out  on  a  data  file. 

A  special  I/O  routine,  TDRR,  Is  used  to  handle  the  equations. 
It  writes  large  records  directly  from  core  v/ithout  using 
buffers.  It  also  has  the  capability  of  recovering  the  records 
In  any  order.  In  this  program  a  record  number  Is  assigned  to 
each  record  as  It  Is  written  in  FORSOL  and  then  we  use  the 
record  numbers  In  the  Inverse  order  to  recover  the  records  In 
the  Inverse  order  in  BAKSUB. 


Subroutine  BAKSUB  controls  the  extraction  of  the 
displacements.  it  reads  the  equations  from  the  data  file  in 
the  reverse  order  to  the  v/ay  they  \^ere  v/ritten  and  completes 
the  Gaussian  reduction  solution  scheme.  The  array  LSTP  is 
used  to  tell  v/hat  displacements  from  the  solution  vector  go 
into  which  displacement  array.  At  the  end  the  displacements 
are  placed  In  the  proper  order  and  written  on  the  displacement 
file,  IDSP.  If  this  is  an  asymmetric  loading  run,  they  v/lll 
also  be  written  on  the  accumulation  file,  I  ACM. 

During  the  forward  solution,  tlie  element  loads  are  written  on 
the  load  print  file,  lELOAD,  for  later  use.  The  element 
material  properties  are  written  on  the  material  property  print 
file,  lEMAT.  The  element  stiffness  matrices  are  written  on  a 
scratch  file,  ISCR3,  and  passed  to  the  stress  module  for  use 
In  the  accuracy  check,  energy,  and  reaction  force 
calculations.  The  stiffness  matrix  and  force  vector  are  also 
written  on  file  ISTFF  for  use  In  the  stress  calculations. 


Iil£  SdELStl  MODULE 


This  module  calculates  stresses  and  strains  for  each  element. 
Calculation  of  energy,  reaction  forces,  and  accuracy  check  Is 
also  done  liere.  If  v;e  have  an  Incremental  loading  run, 
subroutine  IMCLDC  is  entered  to  control  that  function. 

Subroutine  STSMOD  is  the  module  driver;  It  calculates 
dimensions  for  the  basic  arrays  and  passes  control  to 
subroutine  STRDRV. 

Subroutine  STRDRV  calls  the  particular  routines  needed  to 
calculate  the  stresses  and  strains  and  then  calls  the  routines 
for  energy,  reaction  forces,  and  accuracy  check  calculations. 


Ill£  PRINT  MOmi 


This  module  prints  the  data  produced  by  the  other  modules. 
Various  routines  are  called  to  print  the  specific  sets  of 
data . 

The  data  contained  on  files  lELOAD,  lEtlAT,  IDSP  and  I PRMT  are 
all  printed  as  requested.  The  element  loads  are  on  lELOAD. 
The  element  material  properties  are  on  lEMAT.  The  nodal 
displacements  are  on  IDSP.  Stresses  and  strains,  energy, 
accuracy  check,  and  reaction  forces  are  on  I PRMT. 
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Common  SPRINT  Is  used  by  several  of  the  routines  as  a  v/ork 
area.  Its  only  function  is  to  order  the  variables. 


lili  PLOT-ACCUMULATION 


This  module  has  two  functions.  First#  It  produces  deformed 
geometry  plots  on  command  from  the  Input.  Second#  It 
accumulates  and  prints  the  data  from  the  asymmetric  loading 
run.  These  tv/o  functions  may  be  Interspersed  so  that  plots 
can  be  made  of  the  accumulated  data. 

All  plots  are  made  from  data  on  scratch  file  ISCRl.  If  this 
file  has  not  been  generated  by  the  accumulation  function#  it 
must  be  created  by  subroutine  CISCR#  v/hich  is  there  for  that 
purpose. 

Two  sets  of  flip-flop  files  are  used  in  the  accumulation.  One 
set  is  for  displacements  and  the  other  set  is  for  stresses  and 
strains.  As  each  set  of  data  is  read  in#  the  old  summation  Is 
read#  the  new  data  is  multiplied  by  the  proper  trig  function 
and  added  to  it#  and  then  the  accumulated  data  is  written  out. 
At  the  end  of  each  pass#  the  files  change  roles. 

When  making  a  transverse  cut  through  the  body#  the  points  are 
not  necessarily  picked  up  In  the  desired  order  and  must  be 
sorted.  The  routines  that  do  this  are  so  written  that  they 
will  make  use  of  all  available  core  to  avoid  unnecessary  I/O 
operations.  When  merging  is  necessary#  the  number  of  required 
merges  Is  v/ritten  out  on  the  audit  file. 


Ui£  STABILITY  MODULE 


This  module  consists  of  only  one  routine.  Due  to  the 
simplicity  of  the  calculations#  the  lack  of  a  requirement  for 
large  amounts  of  core#  and  Its  different  nature  from  the  rest 
of  the  program#  this  routine  simply  prints  as  it  calculates. 
This  allows  the  program  to  run  In  a  minimum  region  when  this 
option  is  used. 
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film  CHARTS 


Following  are  flow  charts  of  the  basic  program  modules. 
Subroutine  CONTRL  is  the  main  control  routine  and  calls  the 
other  module  drivers.  Flow  charts  of  the  module  drivers  and 
main  logic  routines  follow  in  order. 
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SPECIAL  ROUTINES 


There  are  several  assembler  lan^uace  routines  In  the  program 
v/iiich  require  some  special  attention.  These  are  detailed 
below: 


Subroutine  TORR 

The  purpose  of  this  routine  is  to  provide  a  means  to  v/rite 
large  records  on  disk  and  then  retrieve  them  in  any  order. 

The  routine  Is  written  in  assembler  language  and  uses  BSAM 
reads  and  writes  v/ith  note  and  point  to  accomplish  Its  I/O. 

When  the  open  entry  is  called^  a  getmain  is  Issued  to  get  core 
for  a  table.  This  table  will  be  used  to  hold  the  beginning 
address  of  each  logical  record.  The  index  passed  with  each 
read  or  write  is  used  to  reference  this  table  either  to  store 
the  track  address  for  a  write  or  to  retrieve  it  for  a  read. 

Logical  records  may  be  any  size,  and  the  routine  v/ill  break 
them  up  into  blocks  if  necessary.  The  track  address  of  the 
first  block  in  any  one  record  is  all  tiiat  v/ill  be  saved. 

The  maximum  physical  block  size  is  coded  into  the  routine  as 
the  variable  LIM  and  must  always  reflect  the  track  size  of  the 
disk  drive  being  used.  No  buffers  are  used  since  the  routine 
reads  and  v/rites  directly  from  core. 

There  are  three  entries  into  the  routine: 


CALL  TDOPN  ('FTXX',  CODE,  MAX) 

This  entry  must  be  called  to  initialize  activity  on  the  data 
set.  The  XX  in  FTXX  must  be  replaced  witli  the  two  digits  from 
the  DDNAME  for  the  data  set.  Code  must  l)e  1  at  present.  MAX 
is  the  maximum  number  of  logical  records  to  be  v/ritten.  CODE 
and  MAX  are  both  integers. 


CALL  TDWR  ( REC,  K,  L) 

This  entry  will  write  a  logical  record  from  array  REC  for  K 
contiguous  words.  L  is  the  index  number  for  tills  record;  it 
must  be  unique,  betv/een  1  and  MAX,  and  used  to  request  the 
record  when  it  is  to  be  read.  K  and  L  are  Integers,  while  REC 
may  be  of  any  type. 
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CALL  TDRD  (REC,  K,  L) 

This  entry  will  cause  the  record  assigned  Index  number  L  In 
the  v/rlte  command  to  be  read  Into  array  REC  for  K  words. 

I  This  routine  v;lll  handle  only  one  data  set  at  present. 
Subroutine  lAiSlQ 

The  purpose  of  this  routine  Is  to  rapidly  read  and  write  data 
from  scratch  data  sets  In  binary  form.  The  routine  Is  v/rltten 
In  assembler  language  and  uses  bsam  reads  and  writes  to 
accomplish  Its  I/O.  The  records  are  blocked  Into  physical 
records  of  800  bytes  by  the  routine#  and  a  control  word  Is 
attached  to  each  record#  giving  the  length  of  both  the 
preceding  and  the  following  records.  This  control  word  Is 
used  to  check  record  lengths  on  reads  to  verify  the  data#  and 
to  control  the  backspacing  of  logical  records. 

The  maximum  logical  record  allowed  Is  32767  words.  Records 
will  automatically  be  spanned  over  several  physical  records 
whore  necessary. 

A  maximum  of  50  files  with  FT  numbers  1  through  50  may  be 
specified#  and  any  number  may  be  In  use  at  one  time.  A 
getmaln  Is  Issued  for  1744  bytes  for  each  open  file#  and  the 
area  thus  obtained  contains  the  DCB#  DECO  and  buffers  along 
with  the  buffer  control  parameters.  Either  disk  or  tape  may 
be  used. 

The  routine  has  six  entry  points.  In  the  following 
descriptions#  file#  code  and  count  are  Integers  while  the 
arrays  may  be  any  type. 


CALL  TASS  10  (FILE#  CODE) 

This  entry  must  be  called  to  Initiate  activity  on  each  file. 
File  Is  the  unit  number  and  must  be  from  1  to  50.  Code  Is  1 
If  the  file  Is  to  be  used  for  output  or  2  If  1  Is  to  be  opened 
for  Input. 


CALL  TASWR  (FILE#  ARRAYl#  COUNTl# . . . # ARRAYN#  COUNTN) 

This  entry  v/l  1 1  write  as  many  logical  records  as  there  are 
pairs  of  arrays  and  counts.  File  Is  the  unit  number.  Array 
Is  the  area  to  be  written  from.  Count  Is  the  number  of  4-byte 
words  to  be  written.  Count  must  be  no  larger  than  32767. 
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CALL  TASRD  (FILE,  IRTN,  ARRAYl,  COUMTl, . . . , ARRAYM,  COUNTM) 

This  entry  v;i  1 1  read  as  many  logical  records  as  there  are 
pairs  of  arrays  and  counts.  File  Is  the  unit  number.  Array 
Is  tlie  area  to  be  read  Into.  Count  Is  the  number  of  4-byte 
words  to  be  read.  Input  a  -1  In  I RTII  to  have  TASS  10  dump  on 
an  end  of  file.  If  I RTN  I s  0  on  entry.  It  will  be  left  0  on  a 
normal  return  and  will  be  set  to  +1  If  an  end  of  file  Is 
encountered. 


CALL  TASDCK  (FILE,  I  RTN,  ARRAYl,  COUfJTl,  .  .  . ,  ARR  AYfl,  COUMTf!) 

This  entry  will  read  as  many  loyjcal  records  as  there  are 
pairs  of  arrays  and  counts.  The  records  will  be  retrieved  In 
Inverse  order  to  the  v/ay  they  were  written.  This  accomplishes 
tlie  same  tliln^:  as  a  backspace  -  rend  -  backspace  In  FORTRAN. 
A1 1  parameters  are  the  same  as  for  TASRD. 


CALL  TASREW  (FILE,  CODE) 

This  entry  rewinds  the  file  and  positions  It  as  specified  by 
CODE.  If  CODE  Is  1,  the  file  Is  opened  for  output.  If  CODE 
Is  2,  the  file  Is  opened  for  Input  and  the  first  block  Is  read 
and  checked  and  the  second  block  has  a  read  Issued  on  It. 


CALL  TASCLS  (FILE,  CODE) 

This  entry  closes  the  file  as  specified  by  CODE.  If  CODE  Is 
1,  the  file  v/lll  be  rev/ound  and  all  control  and  btjffer  space 
will  be  freed;  TASSIO  must  be  called  for  this  file  before  It 
can  be  used  again.  If  CODE  Is  2,  the  file  Is  completed  but 
left  open  and  positioned  to  the  end  so  that  It  can  be  read 
using  TASIICK.  TASCLS  must  be  called  with  a  CODE  of  2  before 
TASBCR  Is  called. 


Sut?rQutine  getmain 

This  routine  does  a  conditional  getmain  for  a  very  large 
section  of  core.  Then  when  the  system  returns  and  tells  what 
is  available.  It  frees  the  v/hole  thing,  reduces  the  amount 
available  by  a  specified  amount,  and  does  a  second  getmain  for 
the  reduced  region.  This  allows  us  to  get  all  available  core 
In  our  region  and  still  leave  the  system  v/hat  It  needs. 


Subroutine  ,FR£FRM 


This  routine  Is  designed  to  read  cards  In  a  free-fleld  format. 
The  Input  specifications  are  given  in  the  user's  manual. 

Two  things  should  be  mentioned  here  that  are  not  given 
elsewhere. 

First,  the  routine  uses  what  is  referred  to  as  a  control 
string.  There  Is  one  special  control  string  built  into  the 
program  in  the  material  property  Input  routine. 

A  control  string  may  be  used  to  simplify  the  input  of  a 
sequence  of  logical  records.  Each  digit  in  the  control  string 
represents  a  parameter  in  the  logical  record  and  has  the 
following  values: 

0  -  no  Input  for  this  parameter. 

1  -  this  parameter  will  be  input  in  each  record. 

2  -  input  will  be  made  In  the  first  record  and  will  be 
considered  constant  thereafter. 

For  example,  if  a  ten  value  record  is  available  and  the  2nd 
value  is  constant,  the  4th  and  6th  vary  with  each  record,  and 
no  other  values  v/ili  be  input;  a  control  string  of  C0201010000 
may  be  set  up.  The  first  record  would  then  contain  3  values 
and  each  following  only  2. 

Second,  the  routine  begins  scanning  on  a  record  at  the  end  of 
the  previous  record,  not  at  the  beginning  of  a  card.  The 
difference  then  between  title  and  flag  records  Is  that  title 
records  are  treated  as  cards  and  columns  1-72  are  moved 
directly  into  the  flag  area,  v/hlle  all  leading  blanks  are 
stripped  from  a  flag  record.  This  Insures  that  the  first  word 
in  the  flag  area  can  be  compared  to  a  sequence  of  valid  flags. 

Two  FORTRAN  routines  are  used  by  this  routine:  FRERD  to  read 
file  FT05  and  FRERR  to  put  error  messages  on  FT06. 


Subrout  i  ne  DATJjj 

This  routine  is  used  to  get  the  date  as  a  Julian  integer  and 
the  time  in  cent i seconds.  It  does  this  by  issuing  a  TIME 
macro  in  assembler  language. 


<;ubroutlne  Mli 

This  routine  gives  the  calendar  date  as  an  alpha  string. 
Subroutine  TYME 

This  routine  gives  the  time  of  day  in  hours^  minutes  and 
seconds  as  an  alpha  string. 

Subroutine  XFPTRP 

Literally  X-f loat ing-point-trap.  This  routine  traps  certain 
floating  point  errors  and  takes  corrective  action.  V^/hen  an 
underflow  occurs^  no  error  message  is  given;  the  result  is 
simply  set  to  zero  and  execution  continues. 

Vihen  a  division  by  zero  occurs/  the  result  is  set  to  zero  and 
execution  continues.  When  an  overflow  condition  occurs/  the 
program  terminates  with  a  core  dump  and  a  ijsf»r  completion  code 
of  444. 

This  routine  calls  a  FORTRAN  routine  CTNPWT  to  write  its  error 
messages . 

Subroutine  ELPTIM 

This  routine  obtains  the  CPU  and  wait  time  from  the  system. 
Both  times  are  returned  as  Integers  in  hundredths  of  seconds. 


Subroutine  BLKLD 

This  routine  is  used  to  initialize  a  block  of  core  to  a  given 
value.  The  routine  goes  from  one  address  to,  but  not 
1 ncl udi ng.  a  final  address. 


Subroutine  ABQRIA 

This  routine  creates  an  abend  with  a  core  dump  and  user 
completion  code  of  3276  anytime  it  is  called. 


MIA  files 


There  are  several  basic  data  files  v/hose  format  will  be 
defined  in  the  following  paragraphs.  Not  all  files  will  be 
defined  here  since  several  are  scratch  files  or  their  format 
and  usage  are  relatively  obvious  from  the  program  coding. 

Ill£  Geometry  F.lJ.e 

This  has  the  name  IGEOM.  It  contains  geometry  and  flag  data 
row-wise  for  use  throughout  the  program.  It  Is  created  In 
subroutine  PRGEOM. 

The  record  set  for  each  row  of  nodes  Is  as  follows: 

1.  The  J  Index  of  the  row. 

2.  The  minimum  I  Index  for  the  row. 

3.  The  maximum  I  Index  for  the  row. 

4.  All  of  the  R  coordinates  for  the  row. 

5.  All  of  the  Z  coordinates  for  the  row. 

6.  All  of  the  ICODE  data  for  the  row.  For  each  node 
this  consists  of  16  half-words  whose  usage  Is: 

1  Print  flag 

2  Element  type 

3  Material  code 

4  R  boundary  condition  flag 

5  Z  boundary  condition  flag 

6  Not  used  at  present 

7  e  boundary  condition  flag 

8  Sliding  boundary  condition  flag 

9  Not  used  at  present 

10  Not  used  at  present 
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11  Isoparametric  order  of  Integration 
12-16  Not  used  at  present 

The  Load  File 

This  file  has  the  name  I  LOAD.  It  contains  load  and  boundary 
condition  data  row-wise.  It  Is  created  In  subroutine  PRfiEOM. 

The  record  set  for  each  row  of  nodes  Is  as  follows: 

1.  The  J  Index  of  the  row 

2.  The  load  array.  For  each  node  this  consists  of 
38  double  words  as  follows: 

1  Temperature  of  the  element 

2-5  Pressure  on  faces  1-4  In  that  order 

6-9  Shear  on  faces  1-4 

10-13  Radial  traction  (PR)  on  faces  1-4 

14-17  Theta  traction  (PTH)  on  faces  1-4 

18-21  Axial  traction  (PZ)  on  faces  1-4 

22  Radial  body  force 

23  Theta  body  force 

24  Axial  body  force 

25  Not  used  at  present 

26  Radial  boundary  condition  value 

27  Axial  boundary  condition  value 

28  Not  used  at  present 

29  Theta  boundary  condition  value 

30  Sliding  boundary  condition  value 
31-38  Not  used  at  present 


f-Isure  1  shov/s  the  face  desinnatlons  for  pressure  and  shear 
Koads.  The  numbers  In  circles  give  the  face  codes.  Mote  the 
direction  of  increasing  I  and  J, 

If  we  have  an  I  soi^arametric  element  with  extra  nodes  Indicated 
by  the  X's  in  Figure  1,  the  loads  v/ill  be  applied  at  the 
corners  rather  than  on  the  faces.  The  loads  for  node  A  are  in 
tiie  first  and  fourth  load  locations  of  the  data  for  the  linear 
element  A,  The  loads  for  node  B  are  in  the  first  and  second 
load  locations  for  linear  element  B',  Thus  the  location  in 
the  sequence  gives  the  direction  of  the  load^  the  first 
location  being  in  the  direction  normal  to  face  one#  etc. 


jhL  Matfir.lfl.1  File 

This  file  has  the  name  MATFL,  It  contains  the  material 
property  tables.  It  is  written  In  subroutine  HATP  and  read  in 
subroutine  MATRD, 

The  form  of  the  table  Is  as  follov/s  for  each  material: 

1,  A  material  header  record  consisting  of  5  double 

v/ords  as  follows: 
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1  The  material  number  as  a  floating  point  number. 

2  The  number  of  records  (N)  in  this  material  table 
as  a  floating  point  number. 

3-5  A  24-character  material  name. 


2.  The  material  property  detail  record  consisting  of 

15  double  words.  There  will  be  M  of  these  for 
each  material.  Mote  that  this  is  an  orthotropic 
record. 

1-3  The  three  values  of  Young's  modulus.  The 

isotropic  modulus  is  in  location  1. 

4-9  The  six  values  of  Poisson's  ratio.  The  isotropic 

ratio  is  in  location  4. 


10  The  orthotropic  shear  modulus. 


11-13  The  three  values  of  alpha.  The  isotropic  alpha 
is  in  location  11. 


14  Temperature  if  temperature  dependent;  otherwise 
zero. 

15  Strain  if  strain  dependent;  otherwise  zero. 
Displacement  Li. l.C 


This  file  has  the  name  IDSP.  It  contains  the  displacements  as 
written  during  the  backsubsti tutlon. 


The  file  begins  with  a  flag  record  of  the  following  form; 


1. 


2. 


3. 


One  record 
nodes . 


The  number  indicating  displacements. 

The  mode  if  asymmetric  loading. 

A  bit  flag  defining  the  degrees  of  freedom  in  the 
problem. 

of  the  following  form  is  written  for  each  row  of 


1. 


The  J  index  of  the  row. 


2.  The  displacements  for  the  rov/  as  a  series  of 

arrays  U#  V,  ROT  h  II,  any  of  v/hich  may  be  zero 
length.  For  example,  for  the  ax  I  symmetric 
problem,  only  U  and  VJ  actually  exist. 


29 


/A 


Iii£  Prim  File 

This  file  has  the  name  IPRMT.  It  contains  stresses  and 
strains#  energy#  accuracy  check  and  reaction  forces  to  be 
printed  by  the  print  module.  Any  or  all  of  the  groups  nay  be 
present  depending  on  the  requested  options.  Each  group  Is 
laycd  out  as  follows. 

The  stresses  and  strains  are  preceded  by  a  flag  record: 

1.  The  number  2#  Indicating  stresses. 

2.  The  mode  number  If  asymmetric. 

3.  A  bit  flag  giving  the  degrees  of  freedom. 

The  flag  Is  followed  by  one  record  of  the  following  form  for 
each  set  of  stresses  and  strains  to  be  printed: 

1.  The  1  Index 

2.  The  J  Index 

3.  A  record  containing  the  values.  This  varies  In 
length  and  content  from  element  to  element.  See 
the  respective  stress  or  print  routines  for  the 
detail  for  a  given  element. 

The  energy  Is  preceded  by  a  flag  record  similar  to  the  one  for 
stresses  except  that  the  flag  for  energy  Is  the  number  3.  The 
flag  Is  followed  by  one  record  of  the  form: 

1.  The  mechanical  energy  term. 

2.  The  temperature  constant. 

3.  The  sum  of  the  tv/o  which  Is  the  total  energy. 

The  accuracy  check  data  is  preceded  by  a  flag  record  similar 
to  the  one  for  stresses  except  that  the  flag  Is  the  number  4. 
The  flag  Is  followed  by  ten  records  of  the  form: 

1.  An  I  index 

2.  A  J  index 

3.  A  component  code 

4.  The  value  of  KU  for  the  component  and  node 


so 


I 


5. 


The  value  of  F 


6.  The  difference  KU-F 


The  reaction  force  data  Is  preceded  by  a  flag  record  similar 
to  the  one  for  stresses  except  that  the  flar.  Is  the  number  5. 
The  flag  Is  followed  by  one  record  for  each  node  with  boundary 
conditions  on  It.  The  last  record  Is  followed  by  a  flag 
record  of  all  hexadecimal  01234567  words  to  mark  the  end.  The 
logical  record  is  of  the  form: 

1.  The  I  Index 

2.  The  J  Index 

3.  Three  double  words  giving  the  reaction  forces  In 
the  coordinate  directions: 

Radial 

Theta 

Axial 


S.er.y.lcfi  Life  FI  le 

This  file  has  the  name  ISLF.  It  Is  written  using  FORTRAN' 
binary  I/O.  It  Is  usually  a  tape  and  contains  data  needed  for 
the  service  life  calculations.  Several  sets  of  data  may  be 
stacked  together  on  the  file.  Each  set  of  data  has  the 
following  form. 

A  flag  record  giving: 

1.  The  date  as  a  Julian  date  in  a  one-word  Integer. 

2.  The  time  of  day  In  cent  I  seconds  as  a  one-v/ord 
Integer. 

3.  The  IS-word  title  of  the  run  producing  the  data. 

This  Is  followed  by  one  record  of  the  following  form  for  each 
element : 

1.  The  I  Index 


2. 


The  J  Index 


t 


•'i 


3. 

4. 

5. 

6. 


The  thermal  load 
The  radial  strain 
The  theta  strain 
The  axial  strain 


7. 

8. 

9. 

10. 

11. 

12. 

13. 

14. 


The 

RZ 

shearing 

strain 

The 

Re 

shear i ng 

strain 

The 

Zd 

shear i ng 

strain 

The 

radial  stress 

The 

theta  stress 

The 

axial  stress 

The 

RZ 

shear i ng 

stress 

The 

Rd 

shear i ng 

stress 

15.  The  ZO  shearing  stress 

The  last  record  is  followed  by  a  flag  record  which  indicates 
the  end  of  that  set.  It  Is  composed  of  28  v/ords  of  a 
hexidccinai  01234567. 
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i 

I 

i 
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COMMONS 


Common  /$BEAR/ 


Usage :  This  common  stores 
the  subroutine  BGRID. 

the  bearing  geometry  to  be  passed  to 

Modules:  1 nmod 

Parameters 

name 

Type 

IMIN 

1 

Minimum  1  value  in  the  bearing. 

IMAX 

1 

Maximum  1  value  in  the  bearing. 

JMIN 

1 

Minimum  J  value  in  the  bearing. 

UMAX 

1 

Maximum  J  value  in  the  bearing. 

IVEC 

1 

Vector  flag 

mil 

1 

Dummy  for  word  alignment 

va:jg 

R*8 

Vector  angle  (degrees) 

ITF 

1 

Mot  used  at  present 

1  EQU 

1 

If  this  is  nonzero#  the  nodes 
will  be  equally  spaced  across 
the  bearing  in  the  1  direction 
Instead  of  having  the  first  and 
last  interval  half  the  others. 

IDl 

1 

Not  used  at  present 

ISC 

1 

♦0  spherical  bearing 

♦1  Conical  bearing 

Rl 

R*8 

Inside  radius  of  bearing 

RF 

R*8 

Outside  radius  of  bearing 

Rl  1 

R*8 

Radial  coordinate  of  a  corner 
node 

RIF 

R*8 

Radial  coordinate  of  a  corner 
node 

RFI 

R*8 

Radial  coordinate  of  a  corner 
node 

RFF 

R*8 

Radial  coordinate  of  a  corner 
node 

Zl  1 

R*8 

Axial  coordinate  of  a  corner 
node 

ZIF 

R*8 

Axial  coordinate  of  a  corner 
node 

ZFI 

R*8 

Axial  coordinate  of  a  corner 
node 

ZFF 

R*8 

Axial  coordinate  of  a  corner 
node 

TS 

R*8 

Thickness  of  shim 

NS 

1 

Number  of  columns  In  one  shim 

iSN 

1 

Number  of  shims 

TE 

R*8 

Thickness  of  elastomer 

NE 

1 

Number  of  columns  in  one 

elastomer 

lEN 

1 

Number  of  elastomers 

ID2(7) 

1 

Not  used  at  present 

MS 

1 

Material  number  for  shim 

ME 

1 

Material  number  for  elastomer 

ID5(19) 

1 

Not  used  at  present 

SillMT(50) 

R*8 

Shim  thickness  table 

ELASKSO) 

R*8 

Elastomer  thickness  table 

MATS (50) 

1 

Shim  material  tables 

MATEC50) 

1 

Elastomer  material  tables 
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COMMON  /$BIND/ 


Usage:  This  common  contains  the  Indices  which  define  the 
region  to  be  filled  In  by  the  subroutine  BEARNC. 

Modules:  I nmod 


Parameters: 

Name 

Type 

Usage 

IB 

1 

Beginning  1  Index 

IE 

1 

Ending  1  Index 

JB 

1 

Beginning  J  Index 

JE 

1 

Ending  J  Index 

COMMON 

/ACCOM/ 

Usage: 

ThI  s 

common 

provides  Information  necessary 

accumulate  and  print  displacements,  stresses  and  strains  for 
an  asymmetrically  loaded  body. 

Modules:  Pltmod 

Parameters: 


Marne 

Type 

Usage 

TORZ 

R*8 

Theta  or  Z  cut  location 

NN 

1 

Maximum  number  of  node  points 
In  radial  direction  on 
transverse  cut. 

MODE 

1 

Mode  number  In  the  Fourier 
series. 

FLAG 

L 

First  time  flag  for  data  read 
In  accumulation. 

IDS 

1 

ID-1 

JDS 

1 

JD-1 

-"•r 


KK 

1 

Accumulation  record  counter 

lOPT 

1 

Option  flag  for  axial 
transverse  cut 

or 

KD 

1 

Actual  number  of  node  points  In 
radial  direction  on  transverse 
cut 

COMMON  /ASYM/ 

Usage:  This  common  stores  data  necessary  for 
loading  option. 

Modules:  Driver,  Solmod,  Stsmod 
Parameters : 

Iy££  Usage 

MODE  I  Mode  number 

LMOO  I  Highest  mode  a1 

problem 


COMMON  /CNTRL/ 

Usage:  The 

general  data  Is 

stored  In  this  common. 

Modules: 

Parameters: 

Inmod,  Driver, 

PItmod,  Prtmod,  Solmod,  Stsmod 

Name 

.’voe 

Usage 

TITLE(18) 

R*4 

Title  for  program 

LRUN(20) 

L 

Run  parameter  flags  (In  same 
order  as  run  parameters  In 

Input  document) 

LSOL 

L 

Logical  control 

LSTRS 

L 

parameters  to 

LPRNT 

L 

sped fy  which 

LPLOT 

L 

module  to  enter 

the  asymmetric 


lowed  In  the 
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LOOPl 


L 

LOOP2  L 

LDUM(5)  L 

IMIN  I 

JMIN  I 

I  MAX  I 

UMAX  I 

NTTPS  I 

NREGM  I 

JSHEL  I 

AXSTRA  R*8 

AXSTRS  R*8 

BTEMP  R*8 

COMMON  /CORMAX/ 

Usage:  This  common  stores 
the  amount  of  core  used  in  a 

Modules:  Inmod,  Pltmod 

Parameters : 

Name  lYJIfi 

KMAX  i 

KJL  i 


Iteration  control  logical 

Iteration  control  logical 

Not  used  at  present 

Mini  mum  i 

Minimum  J 

Maximum  I 

Maximum  J 

Mot  used  at  present 
Not  used  at  present 
Not  used  at  present 
Not  used  at  present 
Not  used  at  present 
Base  temperature 

necessary  information  to  compute 
particular  module. 


USflfie 

The  maximum  number  of  v*/ords 
allocated  from  the  array  A  and 
used  in  one  of  the  subroutines 
of  the  module 

Number  of  words  allocated  from 
the  array  A  and  used  in  the 
driver  of  the  module 


COMMON  /OIMEN/ 

Usage:  This  common  stores  dimensions  necessary  to  run  the 
program. 


^7 


Modules:  Inmod,  Driver#  Pitmod#  Prtmod#  Solmod#  Stsmod 
Parameters : 


Name 

T  voe 

Usage 

ID 

1 

The  1  dimension  of  the  problem 

JD 

1 

The  J  dimension  of  the  problem 

MD 

1 

Number  of  rows  of  geometry 
necessary  to  completely  define 
a  row  of  elements 

MD 

1 

Number  of  single  precision 
words  required  for  material 
property  tables 

LL 

COMMON 

L 

/DSPRC/ 

l^ords  remaining  in  array  A 
after  words  have  been  allocated 
to  v/orklng  routines 

Usage: 

Th I s  common 

provides 

work  space  for  those  routines 

which  accumulate  or  print  displacements  for  an  asymmetrically 
loaded  body. 

Modules:  PItmod 

Parameters: 

Name 

Type 

USflfig 

TH 

R*4 

Theta  or  Z  coordinate 

R 

R*4 

Radial  coordinate 

OR 

R*4 

R  displacement 

DT 

R*4 

Theta  displacement 

DZ 

COMMON 

R*4 

/ELPTR/ 

Z  displacement 

Usage: 

This  common 

stores 

pointers#  etc#  to  assemble  the 

element 

stiffness  matrix  into 

the  master  stiffness  matrix. 

TO 


Modules:  Sol mod 


Parameters: 

Name  Tvoe 

NU(3,3)  I 

NW(3,3)  I 

NROT(3,3)  I 

NV(3,3)  I 

MH(3,3)  I 

INDOF(3,3)  I 

MTXS I Z  I 

LIM  I 


Usage 

U  pointers 

W  pointers 

Rotation  pointers 

V  pointers 

H  pointers 

Degrees  of  freedom 

Element  matrix  dimension 

Number  of  degrees  of  freedom 
eliminated  from  element  matrix 


COMMON  /FILES/ 

Usage:  The  common  contains  all  file  designations  except  that 
for  the  TDRR  file#  which  uses  FTOlFOOl#  and  that  for  the 
service  life  file#  which  has  its  own  common.  All  I/O  is  to 
use  the  symbolic  names  as  given  here. 


lliUQg 
I  OUT 
I  PUN 
lAUD 
INPL 

I  ME  IS 
MATFL 
IGEOM 
I  LOAD 


LU  Number  Usage 

6  Printer 

7  Card  punch 

8  Printer  -  audit  file 

12  Plot  and  accumulation  control 
records 

13  Not  used  at  present 

14  Material  properties 

10  Geometry  file 

11  Load  file 
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but 


IDSP 

18 

Displacement  file 

1  PRNT 

19 

Print  file  for  al 1  but 

dl splacements 

ISTFF 

20 

Stiffness  data  passed  from 

stiffness  to  stress  routines 

1  EMAT 

21 

Element  material  properties 

ISCRl 

22 

Scratch  file 

ISCR2 

23 

Scratch  file 

ISCR3 

24 

Scratch  file 

ISCR4 

30 

Scratch  file 

ISCR5 

31 

Scratch  file 

ISCR6/I 

ISDl 

32 

Scratch  file  alro  used  In 

Incremental  loading  option 

ISCR7/ISD2 

33 

Scratch  file  also  used  In 

Incremental  loading  option 

ISCR8 

34 

Scratch  file 

ISCR9 

35 

Scratch  file 

lELOAD 

25 

Element  loads 

1  ACM 

26 

Displacements  and  strains  for 
asymmetric  accumulation 

Note: 

Due  to 

the  fact 

that  some  pairs  of  file  names  are 

Interchanged  In  the  program,  the  names  v/lll  not  always  have 
the  value  given  here.  Also,  different  forms  of  this  common 
will  be  found  In  some  areas  of  the  program;  the  names  may 
change,  but  the  use  of  the  files  Is  as  given  here. 

COMMON  /INCNTL/ 

Usage:  This  common  provides  Information  about  the  correctness 
and  the  type  of  Input  being  read  by  the  subroutine  FREFRM. 


Modules:  Inmod,  Driver 
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Parameters 


Name 

Type 

Usa;te 

IFLG(18) 

1 

Flag  record 

lERR 

1 

Error  return  code 

ICIIT 

1 

Return  code  which  tells  excess 
number  of  values  read 

COMMON  /INCRLD/ 

Usage:  This  common  stores  the  information  necessary  for  the 
control  of  the  Incremental  loading  option. 

Modules:  Inmod^ 

Driver#  Solmod#  Stsmod 

Parameters: 

Marne 

I,^fi 

ILKT 

1 

Counter 

1  LLIM 

1 

Total  number  of  Increments 

1  LPT 

L 

Print  flag 

The  following  8  parameters 
are  applied  Incrementally. 

are  flags  to  indicate  which  loads 

1  LPR 

L 

Pressure 

1  LSil 

L 

Shear 

1  LTllP 

L 

Temperature 

1  LXBF 

L 

R  body  force 

1  LYBF 

L 

Z  body  force 

1  LXfJL 

L 

R  nodal  load 

1  LYML 

L 

Z  nodal  load 

ILND 

L 

Nodal  displacement 

1  LDUM 

L 

Not  used  at  present 

hi 


FACS(IOO) 

R*8 

Incremental  loading  factors 

COMMON  /INCSV/ 

Usage:  This  Is 

an  Incremental 

loading  save  common. 

Modules:  Stsmod 

Parameters : 

Name 

Type 

IG 

1 

File  number  for  IGEOM 

IP 

L 

Print  flag  (LOUT(l)) 

IFLG 

L 

First  time  flag.  Set  TRUE 
the  block  data  and  then 
FALSE  In  IMCLDC. 

COMMON  /ISMAT/ 

Usage;  This  common  Is  used  to  pass  material  properties  for 
the  stress  calculation  on  Isoparametric  elements. 

Modules:  STSMOD 

Parameters! 


Name 

T  voe 

Usar.e 

MU 

R*8 

E/d+NU) 

NU 

R*8 

Poisson's  ratio 

COMMON  /ISOMAT/ 

Usage:  This  common  provides  the  material  counts  and 
numbers  for  Isoparametric  stress  calculations. 

material 

Modulus:  Sol  mod/ 

Stsmod 

Parameters: 

Name 

Type 

Usage 

MTCNT 

1 

Number  of  materials  In 

problem 

MATS(50) 

1 

Material  numbers  In  problem 
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COMMON  /KOMSTF/ 


Usage:  This  common  stores  the  stiffness  matrix  and  force 
vectors  for  the  various  elements  In  the  solution  module. 

Modules:  Solmod 

Parameters : 


Name 

Type 

Usage 

MFULL 

1 

Dimension  of  stiffness  matrix 

LIMNAT 

1 

Number  of  rows  to  be  eliminated 
by  elimination  routine 

S 

R*8 

Stiffness  matrix#  force 
vectors#  etc.#  depending  on  the 
element 

ALFDLT 

R*8 

Thermal  load 

COMMON  /LDCMTL/ 

Usage:  This  common 
the  large  deformation 

stores  the  Information  needed  to  control 
option. 

Modules:  Inmod, 

Driver 

#  Solmod#  Stsmod 

Parameters : 

Name 

.Tyr,e 

Usage 

LDKT 

1 

Iteration  counter 

LDLIM 

1 

Number  of  last  Iteration  to  be 
al lowed 

LDPT 

L 

Flag  for  prInt-out  of 

dl  splacements 

LDCONV 

L 

Flag  for  convergence 

CMVFCI 

R*8 

Convergence  factor 

URF(IOO) 

R*8 

Under  relaxation  factors  for 

Iterations  1-100 


COMMON  /LDKOM/ 


Usage:  This  common  Is  used  to  pass  the  data  required  for  the 
large  deformation  routines. 

Modules:  Soimod,  Stsmod 

Parameters : 


Name 

Type 

Usage 

ALBRll 

R*8 

Non  1 i near 

ALBR22 

R*8 

correct  1  on 

ALBR33 

R*8 

terms 

ALBR13 

R*8 

CC 

R*8 

CBAR 

R*8 

1 

Ef*!! 

R*8 

Strain 

EP33 

R*8 

components 

EP13 

R*8 

EP22 

R*8 

— 

EPl 

R*8 

1 

Principal 

EP2 

R*8 

J 

strains  In  R-Z 

ANGL 

R*8 

Angle  of  principal  stra 

ADIVA 

R*8 

Ratio  of  undeformed 
deformed  area 

fi 

R*8 

iierrman's  variational 

RR(22) 

R*8 

Element  coordinates 

UU(22) 

R*8 

Element  displacements 

iiii 


COMMOri  /LMAT/ 

Usage:  This  common  contains  the  material  properties  for  an 
element . 


Modulus:  Sol  mod 
Parameters : 


Name 

Type 

Usage 

ELMAT(19) 

R*8 

Material  properties  in  i 

orthotropic  form 

ANGL 

R*8 

Orthotropic  angle  of  rotation 

AREA 

R*8 

Area  of  element 

TEMP 

R*8 

Temperature  of  the  element 

COMMON  /MAT PR/ 

i 

i 

Usage:  This 
the  parameters 

common  provides  work  space  for  the  ordering  of  '' 

(material  properties)  within  the  subroutine. 

Modules:  Inmod 

] 

Parameters: 

Marne 

Type 

Usage 

XMAT 

R*8 

Material  number  , 

E(3) 

R*8 

Young's  moduli  j 

XNU(6) 

R*8 

Poisson's  ratios  \ 

G 

R*8 

\ 

Shear  inodulus  ^ 

ALP(3) 

R*8 

Alpha  (these  are  the  thermal  ^ 

loading  coefficents)  ! 

TEMP 

R*8 

Temperature 

r 

STRN 

R*8 

Strain 

ANM(3) 

R*8 

Material  name 
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COMMON  /NRGEE/ 


Usage: 

This 

common 

contains  total  energy  and  variables  used 

for  accuracy 

check. 

Modules:  Stsmod 

Parameters : 

Name 

Type 

Usafie 

TCOM 

R*8 

Element  temperature  constant 

ENTOT 

R*8 

Total  energy 

TOTCON 

R*8 

Total  temperature  constant 

lAC(lO) 

1 

1  indices  of  elements  with  the 
greatest  difference  (KU-F) 

JAC(IO) 

1 

J  Indices  of  elements  with  the 
greatest  difference  (KU-F) 

KAC(IO) 

1 

Flags  which  indicate  the  degree 
of  freedom  for  which  the 
largest  difference  occurred 

ABSR(IO) 

R*8 

KU-F 

CAYU(IO) 

R*8 

KU 

EF(IO) 

R*8 

F 

COMMON  /OUTPUT/ 

Usage:  This  common  contains  the  flags  for  output  options  with 
respect  to  the  entire  output  set. 

Modules:  Inmod,  Driver,  Pltmod,  Prtmod,  Stsmod 


Parameters: 

Name 

Type 

Usage 

LOUT(20) 

L 

These  appear  In 

the 

same 

order 

and  have  the 

same 

meaning  as 

the  parameters 
document. 

I  n 

the 

input 

COMMON  /PAGE/ 


Usage:  The  data  stored  in  this  common  controls  the  printing 
of  page  headings. 

Modules:  Inmod,  Driver,  Pltmod,  Prtmod,  Solmod 
Parameters: 


Name 

NPGS 

I  LNS 
I  PAGE 
LNKNT 


Type 

I 

I 

I 

I 


COMMON  /PLTCOM/ 


Pages  of  output  required  for 
this  print 

Lines  left  on  page 

Current  page  number 

Line  counter 


Usage:  This  common  stores  data  necessary  for  plotting. 
Modules:  Inmod.  Driver.  Pltmod 


Parameters : 

Name 

Type 

IMN 

1 

JMN 

1 

IMX 

1 

JMX 

1 

DSPM 

R*8 

XSCL 

R*8 

YSCL 

R*8 

XM 1  ri 

R*8 

YMIN 

R*8 

Indices  which  define 
the  region  over 
which  we  are 
plotting 

Displacement  multiplier 
X-scale  factor 
Y-scale  factor 
X-mi n imum 
Y-mi n imum 
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XMAX 

R*8 

X-maxlmum 

YMAX 

R*8 

Y-maxInun 

XPS 

R*8 

Paper  size,  X-dl recti  on 

YPS 

R*8 

Paper  size,  Y-dlrectlon 

XTL 

R*8 

Coordinates  for  title 

YTL 

R*8 

location 

ITL 

R*8 

Title  size 

IXY 

R*8 

Orientation  factor 

PTITL(6) 

R*8 

Plot  title 

COMMON  /PRMC/ 

Usage:  The  first  parameter  In  this  common  Is  a  flag  to 
determine  v/hich  of  the  subroutines  In  Prtmod  are  entered. 

Modules:  Prtmod 

Parameters: 


Name  Type 

NFLG(l)  I 


NFLG(2)  or  Mode  I 
NFLG(3)  I 


Usar.£ 

Flag 

1- dl splacements 

2- stress 
5-energy 

4- accuracy  check 

5- reactlon  forces 

Mode  number  If  asymmetric 

Bit  flag  giving  the  degrees  of 
freedom  In  the  problem 


48 


COMMON  /ROWEL/ 


Usage:  This  common  is  used  in  the  solution  module  to  pass 
indices. 


Modules:  Solmod# 

Parameters: 

Stsmod 

Name 

lYJlg 

UsaAfi 

IMN(3) 

1 

Minimum  1  index  for  the  rows  of 
nodes  in  core 

IMX(3) 

1 

Maximum  1  index  for  the  rows  of 
nodes  in  core 

1 

1 

1  index 

J 

1 

J  index  for  a  row  of  elements 
in  core 

J6 

1 

J  index  of  row  of  elements  in 
core 

LTYPE 

1 

Element  type 

IWD 

1 

Number  of  nodes  across  an 

element 

ji/n 

1 

Number  of  rows  of  nodes 
required  for  the  J6th  row  of 
elements 

12 

1 

Same  ns  IWD  for  this  version 

J2 

1 

Same  as  JWD  for  this  version 

MUMPAR 

1 

Maximum  number  of  degrees  of 
freedom  in  problem 

COMMON  /SLFILE/ 


Usage:  This  common  contains  the  designation  for  the  tape 
storing  the  information  needed  in  the  service  life  routine. 

Modules:  Driver#  Pltmod#  Prtnod 


Parameters: 


flame 

Iy£& 

Usar.e 

ISLF 

1 

Service  life  file.  Set  to  9 

for  FT09F001  In  the  block  data 

COMMOU  /SOLDIM/ 

Usage;  This  common  stores 
solution  module. 

the  Important  dimensions  for  the 

Modules:  Solmod 

*  i 

Parameters ; 

}- 

i 

Name 

lyp.c 

Usage 

• 

i 

K4 

1 

Maximum  number  of  equations  to 
be  reduced  at  one  time 

( 

K5 

1 

Semi -bandwl dth 

1 

1 

I 

K6 

1 

Number  of  equations  for  a  row 
of  elements 

1 

i 

t 

IH 

1 

Number  of  equations  to  be 
reduced  at  one  time  for  a  row 
of  elements 

MF 

1 

Column  limiter  for  reduction 

1 

NCOL 

1 

Semi -bandwl dth  for  a  row  of 
elements 

h 

i 

NS 

1 

Same  as  K6  at  present 

1 

i 

NREC 

1 

Record  count  on  the  solution 
file 

i 

IL 

1 

Minimum  1  for  a  row  of  elements 

i 

i 

1 

lU 

1 

Maximum  1  for  a  row  of  elements 

! 

K3 

1 

Required  size  of  solution 

vector 

1 

K9 

1 

Number  of  degrees  of  freedom  on 
a  row  of  nodes 

1  1 

1 

NDOF 

1 

Number  of  degrees  of  freedom  on 
a  node 

1 
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COMMON  /SPRINT/ 


Usage;  This  common  provides  work  space  for  the  ordering  of 
the  parameters  (stresses  and  strains)  within  the  subroutine. 

Module;  Prtmod 

See  the  using  routines  for  structure. 

COMMON  /ST  I  MOD/ 


Usage;  This  is  the  common  for  the  stability  module. 
Modules;  Inmod^  Stamod 
Parameters ; 


Marne 

X 

Y 

A 

B 

Z 

C 

COMVF 

MAXITR 


Usage 

These  variables 
are  the  Input  to 
the  stability 
mod  u 1 e 


ID  I  Dummy  for  boundary  alignment 

XNOS  R*8  Number  of  shims 

OSF  R*8  Old  shape  factor 

COMMON  /STOUT/ 

Usage;  This  is  the  stress  output  common. 

Modules;  Stsmod 
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Parameters: 


Name 

LEN 

INK 

ENRG 


I 

I 

R*8 


COMMON  /STS DIM/ 


Usage 

The  number  of  values  becinning 
with  ENGR  to  be  written  on 
I  PR  NT 

The  Index  Into  the  array  where 
the  Incremental  loading  accumu¬ 
lation  begins 

Element  energy 

The  remainder  of  this  common 
will  vary  In  length  and  content 
according  to  the  stress  routine 
Involved. 


Usage:  This  common  stores 
the  stress  module. 

Modules:  Stsmod 

Parameters : 


Name  Tvoe 

NDOF  I 

I  DOF  I 

KS  I 

IPTR(60)  I 


an  array  of  parameters  needed  by 


Usage 

Number  of  degrees  of  freedom 

A  bit  flag  to  Indicate  the 
degrees  of  freedom  In  the 
problem 

Number  of  displacements  on  a 
row  of  nodes 

An  array  of  pointers  used  to 
give  the  relation  between  the 
element  stiffness  matrix  and 
the  assembled  stiffness  matrix 


COMMON  /STSMAT/ 

Usage:  This  comnran  contains  the  material  properties  for  an 
element . 
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Modules:  Stsmod 

Parameters : 

Name 

T.Y.Pe 

Usar,g 

ELMAT(19) 

R*8 

The  material  properties  In 

orthotropic  form 

ANG 

R*8 

Orthotropic  angle  of  rotation 

AREA 

R*8 

Area  of  element 

TEMP 

R*8 

Temperature 

COMMOM  /STSSTF/ 

Usage:  This  common  stores  the  stiffness  matrix  and  force 

vectors  for  the 

various 

elements. 

Modules:  Stsmod 

Parameters: 

Marne 

Tvne 

Usaite 

MFULL 

1 

Dimension  of  stiffness  matrix 

L 1 MNAT 

1 

number  of  rows  to  be  eliminated 
by  the  elimination  routine 

3(20,24) 

R*8 

The  stiffness  matrix,  force 

vectors,  etc.,  depend  in, on  tlie 
element 


COMflOfl  /SRTCOM/ 

Usage:  This  common  stores  Information  relative  to  the  sorting 
of  records  and  the  merging  of  files. 

Modules:  Pltmod 

Parameters : 

flams  usa^.e 

L5  I  Block  length  vs  records  for 

input 
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N5 

1 

M6 

i 

NM 

1 

m 

1 

LEN 

1 

MAX 

1 

COMMON 

/STRNRC/ 

Usage : 

This  common 

prov i des 

v/h  i  ch 

accumulate 

or 

print 

asymmetrically  loaded 

body. 

Modules 

:  Pltmod 

Parameters : 


Name 

Type 

TH 

R*4 

R 

R*4 

MU 

R*4 

NU 

R*4 

EPR 

R*4 

EPT 

R*4 

EPZ 

R*4 

GAMRZ 

R*4 

GAMRT 

R*4 

GAMZT 

R*4 

H 

R*4 

Number  of  records  per  file  on 
files  ISCR5  and  I SCR6 
respect ively 

Number  of  ner/jes  required 

Maximum  number  of  records  that 
can  be  held  in  core 

Humber  of  words  per  record 

Total  number  of  records 

work  space  for  those  routines 
stresses  or  strains  for  an 


UsafiQ 

Theta  or  Z  coordinate 
Radial  coordinate 
Shear  modulus 
Poisson's  ratio 
R-stra i n 
Theta  strain 
Z-stra i n 
RZ  shear 
R-theta  shear 
Z-theta  shear 
Merrman's  variational 


TiJERM 


Thermal  load 


R*4 

COfinOH/TROUBL/ 

Usage:  This  common  contains  an  error  flag  for  termination  of 
program. 

Modules:  Inmod,  Driver,  PI tmod,  Prtmod,  Solmod,  Stsmod 
Parameters : 

Marne  Type  IJsaire 

IRTRM  I  Error  flag 

«  0  NORMAL  RETURN 
.GT.O  TERMINAL  ERROR 
.LT.O  TRY  NEXT  CASE 

COMMOil  /WORK/ 

Usage:  This  is  a  work  common  which  also  stores  material 
properties  and  loads  for  an  asymmetric  loading  element. 

Modules:  Sol  mod 

Parameters : 


Marne 

Type 

Usaite 

EE 

R*8 

Young's  modulus 

PNU 

R*8 

Poisson's  ratio 

BFRL 

R*8 

BFTHL 

R*8 

Body  forces 

BFZL 

R*8 

THRL 

R*8 

Thermal  load 

P 

R*8 

Pressure 

TAU 

R*8 

Shear 
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PTH 

R*8 

PZ 

R*8 

RS(3) 

R*8 

ZS(3) 

R*8 

AR 

R*8 

t 

4 


Traction  loads 

Triangular  element 
coordinates 

A:  ea  of  a  triangular  element 


£&RQ£  CODES 


There  are  several  user  completion  codes  that  may  be  given  by 

the  program.  These  are  given  here  as  an  aid  in  debugging. 

100  Issued  by  TDRR.  An  end  of  file  was  encountered  on 

a  read. 

200  Issued  by  TDRR.  A  residual  count  error  was 

detected  on  a  read. 

444  Issued  by  XFPTRP.  An  exponent  overflow  condition 

was  detected. 

500  Issued  by  GTMAIM.  Insufficient  core  available  to 

continue  the  run. 

1000  Issued  by  TASSIO.  A  permanent  I/O  error  has 

occurred  on  a  data  set. 

1010  Issued  by  TASSIO.  An  error  has  been  detected  in 

the  TASSCB.  This  may  be  caused  by  a  bad  file 
number  or  by  the  TASSCB  being  destroyed, 

1020  Issued  by  TASSIO.  Attempted  to  open  a  file  that 

v/as  already  open. 

1030  Issued  by  TASSIO.  Error  In  a  getmain  required  to 

open  a  file. 

1040  Issued  by  TASSIO,  An  error  occurred  in  the  open 

routine, 

1110  Issued  by  TASSIO.  An  attempt  was  made  to  v/rlte  N 

words  where  N  is  less  than  one. 

1120  Issued  by  TASSIO.  Attempted  to  v/rlte  on  a  file 

that  was  not  open. 

1200  Issued  by  TASSIO.  An  end-of-file  v/as  encountered 

and  IRTM  v/as  set  to  -1. 

1210  Issued  by  TASSIO.  Attempted  to  read  M  words  where 

M  was  less  than  zero. 

1220  Issued  by  TASSIO.  Attempted  to  read  a  file  that 

was  not  open. 


1230 

1310 

1320 

3276 


Issued  by  TASSIO.  VIrong  length  record.  The  number 
of  words  requested  did  not  match  the  number  of 
words  In  the  record. 

Issued  by  TASSIO.  Tried  to  do  a  backread  on  a  file 
that  was  not  open. 

Issued  by  TASSIO.  Wrong  length  record  on  a 
backread.  The  number  of  words  requested  did  not 
match  the  number  of  words  In  the  record. 


Issued  by  ABORTA.  This  routine 
various  places  In  core  where  an 
requiring  termination  with  a 
detected.  The  traceback  will  show 
occurred. 


Is  called  at 
error  condition 
dump  has  been 
where  the  error 
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aViHLA-Y 


Only  the  philosophy  of  the  overlay  structure  v/ill  be  given 
here.  The  detail  of  the  structure  can  be  readily  seen  from 
the  linkedit  list  or  a  list  of  the  overlay  cards. 

We  have  attempted  to  set  up  the  overlay  In  a  modular  form  to 
match  the  structure  of  the  program.  This  v/ay  only  the  modules 
that  are  actually  used  will  be  loaded. 

Within  each  module  we  have  tried  to  group  commons  and  routines 
together  so  that  only  the  needed  routines  will  be  loaded  and 
also  so  that  the  chains  are  as  nearly  equal  in  length  as 
possible  under  the  other  restrictions. 

This  has  resulted  in  a  highly  overlayed  program  v/I  th  a  fairly 
flat  structure.  The  solution  module  Is  the  longest  module# 
and  care  should  be  taken  that  the  program  Is  not  unduly 
lengthened  by  modifications  to  this  module. 
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